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Abstract 

The proton size, specifically its charge radius, was thought known to about 1% accuracy. Now 
a new method probing the proton with muons instead of electrons finds a radius about 4% smaller, 
and to boot gives an uncertainty limit of about 0.1%. We review the different measurements, some 
of the calculations that underlie them, some of the suggestions that have been made to resolve the 
conflict, and give a brief overview new related experimental initiatives. At present, however, the 
resolution to the problem remains unknown. 
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1 Introduction 

The proton radius conflict is a prime problem in proton structure physics. To state the problem is 
simple: we measure the proton radius using electrons, and we measure the proton radius using muons, 
and we get incompatibly different answers. 

Measurements made using electrons come either from electron-proton scattering or from atomic 
spectroscopy, where there are small but measurable shifts in the hydrogen spectrum due to the proton 
size. The proton charge radius obtained from different electron measurements substantially agree, and 
gave a result with an uncertainly of order 0.6%. 

On the muonic side, there is only one set of measurements, but it is one with a very small uncertainty 
limit. The announcement came in 2010 (^, with follow-up in 2013 [^, of the proton size measured from 
its effect on the muonic hydrogen spectrum, specihcally the 2S-2P Lamb shift. The muon, because of 
its mass, orbits closer to the proton than an electron, and it energy levels are more strongly affected by 
the proton size. This allows a proton radius measurement with an error bar more than 10 times smaller 
than obtained from the collected electronic measurements, despite the limited lifetime of the muon. 

The proton radius from the muonic Lamb shift is a startling 4%, or 7 of the current electronic 
standard deviations, smaller than the previously accepted result. 

The reasons for this are not yet clear. We will make some comments here about what has been 
thought about, beginning with a brief description of the electron based experiments followed by some 
remarks about how the muonic Lamb measurements are carried out to give the charge radius. We will 
then comment on some (but far from all) of the theory calculations needed to interpret the experiments. 
Much work has been done checking the corrections involved in pulling out the charge radius from the 
muonic experiment, but nothing big enough to affect the result that has been found. Then we will 
comment at some length about suggestions that have been made to reconcile the two measurements, 
including speculation about anomalous QCD corrections, speculation about connections to physics 
beyond the standard model, and questions raised about the interpretation of the experiments using 
electrons. The importance of the problem has led to a number of experiments to check aspects of the 
conflict, and we will catalog some of them in anticipation of their success. We will offer some concluding 
remarks, but the main concluding remark, unfortunately or otherwise, is that the problem remains. 

Previous reviews are found in [^|^. 


2 The electronic measurements 

Measurements of the proton radius using electrons, either by electron-proton scattering or by careful 
measurements of atomic energy levels, did not originate the proton radius puzzle. The various electron 
based measurements of the radius are in agreement and it was not until the muon based measurement 
appeared that there was a problem. However, we will start with a discussion of the electron based 
measurements in order to show our notation and delineate one side of the conundrum. 
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We start with electron scattering. For understanding some definitions, it is even useful to start 
non-relativistically, where scattering electrons off pointlike protons has the Rutherford cross section 


da 

dVt 


point 


a 


AE sin\e/2) 


( 1 ) 


in the lab system, where E is the energy of the incoming electron, 9 is its scattering angle, a is the fine 
structure constant, e^/(47r). The nonrelativistic cross section result off an extended target is modified, 
but just becomes 
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where Q is the momentum transfer, the difference between the incoming and outgoing projectile mo¬ 
menta, and G{Q^) is the “form factor”, given nonrelativistically by 


G(Q2) W j d\ IV’WP, 


(3) 


where 'ip{r) is a wave function describing the density of the target. It is easy to expand the form factor 
at low and obtain, 

G(Q") = 1 - hr^)Q" + ..., (4) 

0 

where (r^) is the rms radius squared. 

For a electron-proton scattering at the energies under consideration, relativity matters. The lowest 
order calculation is one photon exchange and proceeds from the Feynman diagram in Fig.|^ The proton 
vertex involves the matrix element of the electromagnetic current between incoming and outgoing states 
of definite momentum, 

= h(p') u{j )), (5) 

where a^y = (*/2)[7^, 7 i,], q = p' — p, E^ and Ep are the Dirac and Pauli form factors, spin degrees of 
freedom are tacit, and = —{p' — pY- 

The electric and magnetic form factors are linear combinations of Ep and Fp, 

Gm = Fp + Ep , (6) 

Ge = Ep-tEp, ( 7 ) 

and the differential cross section to leading order in perturbation theory is given by 

X jj^(G|(7) + Ti,(Q^)). (8) 

AM^ ’ 

n 

= 1 -h 2(1 -h r) tan^ - , (9) 

and and da/dfl\j^g is the Mott or “no-structure” cross section (specifically, the hypothetical cross section 
for scattering an unpolarized spin-1/2 electron from a pointlike spin-0 target), 
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Here and above. 
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Figure 1: Spectrum of ordinary hydrogen (not to scale). 


where E and E' are the incoming and outgoing electron laboratory energies. Obtaining the form factors 
by measuring the differential cross section at a variety of and scattering angles is termed a Rosenbluth 
separation. 

Relativistically, Eq. ([^ is an approximation valid if the target’s Compton wavelength is small com¬ 
pared to its size. For a proton the latter statement is marginal. Instead we dehne the form factors from 
the matrix element given above, and for the electric and magnetic radii, we promote the former derived 
result to a dehnition. 
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( 11 ) 


with a similar relation between Rm and the magnetic form factor. 

The most accurate experiments for obtaining the form factors at low and moderate are done at 
Mainz and reported in [^. They quote their results by htting form factors using a variety of different 
htting functions to their measured differential cross sections. Their hts lead to 


= 0.879(8) fm, (12) 

where, following others [^, several uncertainties are combined into a single uncertainty limit. 

The other way to obtain the proton radius using electrons is to measure the energy levels, the Lamb 
shift but not only the Lamb shift, in ordinary electronic hydrogen. The proton radius measurements 
made using ordinary hydrogen are quite remarkable because the proton radius effects are very small. 
Energy level splittings can be measured so precisely, and the proton size independent energies can be 
calculated so reliably, that proton size dependent terms can be isolated. 

As a reminder, the proton energy spectrum is illustrated in Fig. The Figure is not to scale. (That 
the Lamb shift—the splitting between the 2Si/2 and 2Pi/2 levels—is about 10% of the 2P hne structure 
splitting is the only relative splitting that is about right.) 

The energy levels, at least the S-state levels, are affected by the hnite proton size. The size dependent 
energy shift was worked out non-relativistically in [^, and is to leading order in perturbation theory 
(in modern notation) 

AEi = —102(0) |R|, (13) 

where 0(0) is the electron wave function at the origin in coordinate space, 0^(0) = (m^a)^/(vrn^), where 
nir is the reduced mass and n is the principal quantum number. Interestingly, the result depends only 
on the RMS proton radius squared and not on the detailed distribution of matter within the proton. 
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Relativistic considerations may be examined in , and the resnlt is that the formnla is the same with 
the R\ appearing being precisely the R\ dehned from the derivative of the electric form factor. 

For the discussion, the measurements of the charge radius from atomic energy splittings can be 
divided into small splitting measurements, which are measurements of the corrections to the Lamb 
shift, all levels having the same principal quantum number, and big splitting measurements, which are 
corrections to splittings between levels with different principal quantum numbers. Every splitting has 
a QED term, which is what the splitting would be if the proton were pointlike, and the QED term is 
written as some factor times the Rydberg constant, which is the leading order binding energy of the 
ground state hydrogen atom in a inhnite proton mass limit (and in some units is a^meC^/2). Also, every 
splitting has a proton radius term, which is noticeable only if a S-state is involved. Generically, 

AE = a ■ Ryd + bR^ . (14) 


The coefficients a and b are known; for a, which must be known to high accuracy, this entails extensive 
calculation. 

For the small splitting measurements, the Rydberg constant is accurately enough known from other 
sources to subtract the Rydberg constant term from AE and obtain the R\ term. For the big splitting 
measurements, the R\ term is tiny compared to the main term and the Rydberg constant is not 
independently known well enough to accurately obtain the R% term from a single measurement. The 
solution is to measure two splittings, and solve the elementary pair of simultaneous equations to obtain 
both the Rydberg constant and Re- To repeat: the high precision value of the Rydberg constant is 
obtained from the very same atomic energy level experiments that measure the proton radius. One gets 
the Rydberg constant to about 5 parts per trillion (10^^) and the proton radius to a bit better than 1% 
when all averaging is done 
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A summary of the experimental hydrogen atom results in graphical form is shown in Fig. The 
top three results are Lamb shift measurements. The next 12 are from pairs of larger splittings. The 
1S'-2S' splitting is very accurately measured and is always used as one member of the pair; after that 
there is wide selection of level splittings. In general the spectral lines have a natural width because the 
states can decay, and to obtain the quoted accuracy, the centers of the lines must be located to a part 
in 100 to a part in 1000 of the natural width. One notices immediately that none of the results is a 1% 
measurement. The sub-1% uncertainty limit comes from averaging many experiments. 

The proton radius measurements using electrons are summarized (as part of a extensive summary of 
physical and technical data) by the Committee on Data for Science and Technology (CODATA), which 
has been updating results every 4 years. When the first muonic results appeared, the 2006 CODATA 
results were available. These included some electron scattering results with about 2% uncertainty 
limits 11 , but were dominated by the atomic measurements and reported Re = 0.8768(69) fm. The 


up-to-date number is the CODATA 2010 value, which by the CODATA protocol is based on all data 
available through the end of 2010 and so was not available in 2010 itself. The current CODATA average 
using only electronic spectroscopy data is 0.8758(77) fm 10 . With the Mainz electron scattering result 
included, CODATA finds the overall result (Tol, 


Re = 0.8775(51) fm. 


( 16 ) 


The electron based proton radius results are mutually compatible. Within the available accuracy, 
the electron based experiments by themselves create no problems. However, there came the tempting 
prospect of increasing the accuracy by a factor of ten, and pursuing that prospect led to a great surprise 
and the present day conundrum. 
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Figure 2: Proton radius from ordinary hydrogen level splitting measurements. 


3 The muonic Hydrogen experiment 


The muon is, as far as we know today, like an electron in all respects except in its mass or things directly 
related to the mass, in particular, one can make a hydrogen atom using a negative muon in place of an 
electron. The muon is about 200 times heavier than the electron, and in an atomic bound state it lies 
about 200 times closer to the proton than the electron, so proton size effects are greatly magnified. 

Despite the muon’s limited 2.2 /zs lifetime, it was anticipated that the larger impact of the proton 
size on the energy levels would allow a 0.1% measurement of the proton charge radius. In 2010, the 
promise was more than realized, as a collaboration of experimenters announced their first result [^, 
based on measuring one 2S-2P energy splitting. The same collaboration has more recently published 
measurements of a second splitting, and slightly improved the analysis of the hrst, with the current 
result [^, 

Re = 0.84087(39) fm (16) 


This is 4% smaller than the CODATA radius, and is a 7cr discrepancy. 

The splittings measured in the muonic experiment are splittings in the 2S-2P system, i.e., they are 
Lamb shift measurements, and so are the analogs of the “small splitting” measurements described in the 
electron section. The experiment did not prove easy to make work, but once working, the description 
of it in outline is straightforward. 

The experiment is done at the Paul Scherrer institute in Switzerland by the CREMA (Charge Radius 
Experiment with Muonic Atoms) collaboration and begins with a beam of muons that are slowed and 
then sent into a hydrogen target, where they are hnally stopped and captured on hydrogen atoms in a 
variety of orbits. The muons cascade down very quickly, and nearly all the muons cascade down to the 
15 state, with few stuck in the 2S state, which is metastable. Only about 1% of the muons are in the 
2S state, but these are the ones of interest. After a short delay, to allow the aforementioned cascades to 
hnish, the experimenters shine a tunable laser on the hydrogen, and if the frequency of the laser light 
is just right, there will be a transition to a 2P state; see the left hand panel of Fig[^ 

The 2P state decays very quickly to the ground state, emitting an X-ray, and one has an X-ray 
detector to signal success in tuning the laser correctly. Then the laser frequency gives the desired 
energy splitting. The value of h, incidentally, is only known to a bit better than 7 figures 10 ; however. 
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Figure 3: Energy levels in muonic hydrogen. On the left, the 2S-2P transition is induced by 
a tuned laser, with success indicated by an X-ray from the 2P-1S transition. On the right is 
a detail of the n=2 energy levels of muonic hydrogen. 



Figure 4: Two-photon exchange correction. 


6 hgures suffice for the case at hand. 

The two lines that are measured are indicated in the right hand panel of Fig. where the 25 and 
2P states of muonic hydrogen are shown in more detail. (The subscripts on 2S and 2P indicate the 
summed angular momentum of the orbital motion and electron spin, and F indicates the total angular 
momentum, orbital motion plus electron spin plus proton spin.) The P-states are not measureably 
affected by proton size effects, but the 5-states are, and the proton size affects the hyperhne splitting 
as well as to the overall 5-state energy level. The strict dehnition of the Lamb shift is the 25i/2- 
2 Pi/ 2 splitting with hyperhne splitting effects removed. (The last is easy to do, since the spin-spin 
interaction moves the P = 1 states up by 1/4 of the overall splitting and the P = 0 states down by the 
complementary amount; for more detail see [I^.) 

The two splittings measured are the original 25^^"^^ to 2P^i^‘^ and the more recent 2S^i^^ to 2P^2 ^ 
with splittings in energy units denoted hut and hug, respectively. With two splittings one can obtain 
both the Lamb shift and the hyperhne splitting. The former is the main goal of the experiment, but 
the hyperhne splitting also has proton structure ehects, this time dependent on the Zemach radius 13 


to be discussed below, and the Zemach radius can be determined from this experiment, albeit not with 
the same precision as the charge radius. 

The Lamb shift and hyperhne splitting are 


1 3 

Pl = — 8.8123(2) meV, 

-Fhfs = hug — hut + 3.2480(2) meV, (17) 


where the numerical terms follow from reliable (proton structure independent) calculations of the P-state 
energy splittings. 
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Exp er iment ally, 


= 202.3706(23) meV. 


(18) 


From theory, 


E 


= 206.0336(15) - 5.2275(10)i?| + 


TPE : 


(19) 


where the energies are in meV and Re is in fm. The hrst numerical term includes a number of 
corrections that are independent of the proton charge radius, and the last term is the two-photon 
correction. Fig. The latter will be discussed in its own section below. For now we will mention that 
if the blob between the two photon attachments on the proton side is restricted to being a proton itself, 
and if the calculation is done non-relativistically, the result is proportional to the average separation- 
cubed of the proton charge, denoted can End it so written in older writing {e.g., jl 

Today one uses the value 

= 0.0332(20) meV 


14 


)■ 


( 20 ) 


based largely on calculations in 15 17 and improved in 18 . Equating the theoretical and experimental 


results for the Lamb shift yields the previously stated result, Fq. (16). 

Alternatively, had we used the CODATA 2010 charge radius to calculate the expected Lamb shift, 
there would have been an 320 /reV discrepancy from the experimental value. Either we accept the 
smaller radius or we must hnd about 320 /reV—about 10 times the value of Ftpe— of further corrections. 


In addition, the muonic hydrogen value of the hyperhne splitting, Eq. (17), is 


= 22.8089(51) meV. 


( 21 ) 


This is in good accord with existing calculations of the hyperhne splitting 19 21 . For example. Ref. 21 


gives Fhfs = 22.8146(49) meV. In that work, the Zemach radius is obtained from form factors ht to 
electron scattering experiments. Alternatively, since the proton structure corrections are about 0.7% 
of the hyperhne splitting in muonic hydrogen (compared to about 40 ppm in electronic hydrogen), one 
may ht the Zemach radius from the experimental data, to obtain 

Rz = 1.082(37) fm, (22) 

which this time in line with the electron scattering values, albeit for this quantity there is a 3.4% 


uncertainty 22 . 


3.1 Muonic deuterium 

There being a problem involving ordinary muonic hydrogen with no clear explanation, one may seek 
additional information about the discrepancy by studying other systems. The CRFMA group has also 
measured the Lamb shift in muonic deuterium. As of this writing, the results have been discussed in 
open conferences but not published in a refereed journal. Accordingly, we shall make some comments 
on the expectations and theory involved in this measurement, and only a little about the actual results. 
Three splittings in the muonic deuterium 2S state have been measured, and are diagramed in Fig. 
Since the deuteron is spin-1, the hyperhne splittings triple the states other than the J = 1/2 states, 
and the three measured splittings are from the 2 S'i/ 2 ,f= 3/2 to the 2 P 3 / 2 ,f= 5/2 and from the 2 S'i/ 2 ,f=i /2 to 
both the 2P3/2 ,f= 3/2 and 2P3/2 ,f=i/ 2 states. The measured 251/2, f=i /2 to 2P3/2 ,f= 3/2 splitting is about 
215 meV and the others are close to this value. 

As a simple note, the deuteron is much larger than the proton, so the size ehect is correspondingly 
larger. If the radii had their CODATA 2010 values, the proton hnite size would shift the muonic 25 














Figure 5: A (not to scale) diagram of the n = 2 energy levels of muonic deuterium, with 
indication of the three level spacings that have been measured. The 251/2, f=i /2 to 2P3/2 ,f= 3/2 
energy splitting is about 215 meV. 


state up by 4.0 meV, whereas the deuteron hnite size would shift the same state’s energy up by 27.9 
meV. We already know the discrepancy in the proton case corresponds to lowering the 25 energy by 
about 320 /reV. Without yet understanding the mechanism, we cannot say what the discrepancy would 
be in the deuteron case. For informational purposes, one explanation (there is more discussion in a 
later section) involves a new particle that couples more to muons than to electrons, and also couples to 
protons. If it coupled to the deuteron with the same strength as it coupled to protons (a dark photon 
related explanation would satisfy this criterion), then the energy discrepancy in the deuteron case would 
be the same except for a reduced mass correction, or about 385 /reV. 

Without committing ourselves to this or any other model, we will take 385 peV as a benchmark 
expectation for the energy discrepancy. That is, the normal deuteron size effect will be much larger 
than for the proton, but the discrepancy from the expected value will be about the same as for the 
proton. This benchmark tells us that any uncertainties in expectations or uncertainties in theoretical 
corrections will have to be small compared to 385 peV when translated into energy shifts, if they are to 
be useful. 

A first comment concerns the direct measurement of the deuteron radius in known experiments using 
electrons. One can scatter electrons, obtain a deuteron charge form factor Gc{Q^) and then obtain the 
deuteron charge radius. The result, from just the scattering experiments, is 


Rd = 2.130(10) fm. 


(23) 


or = 4.537(43) fm^. The uncertainty in translates into an uncertainty of about 270/reV in 
the muonic deuterium 25 energy, so the current e-d scattering experiments are unfortunately only of 
marginal use in predicting muonic deuterium energy levels to the required accuracy. New electron- 
deuteron scattering experiments are underway at Mainz to signihcantly improve this situation 23 , and 
we look forward to their completion. In the meanwhile there is another path to take. 

Very accurate measurements of the 15-25 splitting for both ordinary electronic hydrogen and deu¬ 
terium are possible, and the isotope shift leads to very accurate results for the differences of the radii- 
squared. 


Rl-Rl = 3.82007(65) fml 


(24) 
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This combined with the CODATA 2010 value of the proton radius leads to the CODATA 2010 value of 
the deuteron radius, which is 


Rd = 2.1424(21) fm, 


(25) 


where one notes the uncertainty limit is about 5 times smaller than for the radius obtained from 
scattering alone. The resulting uncertainty in the predicted deuteron size dependent energy shift in 
muonic deuterium is now 55 /xeV, which is acceptably small. 

The isotope shift measurement also leads to an alternative estimate the discrepancy from the CO¬ 
DATA radius value that might be expected in the muonic deuterium Lamb shift. If the isotope shift is 
the same for the muonic atoms as for the electronic atoms, the 320 /reV discrepancy in the proton case 
translates into the same discrepancy in the deuteron case, up to reduced mass effects, which is to say 
it leads to the same 385 /reV benchmark discrepancy mentioned above as following from a particular 
model. 

A serious consideration in the muonic deuterium Lamb shift is the theoretical correction from the 
polarizability correction, or equivalently from the two-photon corrections. Fig. In the proton case 
these were important because the theory uncertainty in the calculation of these corrections 15 -18 was 
the largest source of uncertainty in the analysis of the measurement. However, the actual total size of 
the correction was small enough that even a large error in its calculation would not affect that fact of 
a discrepancy between the electronic and muonic radius measurements. The same is far from true in 
the deuteron case. Inelastic contributions to the polarizability requires an energy input of only the 2.24 
MeV deuteron binding energy, whereas a proton requires at least a pion mass of energy for a strong 
interaction inelastic contribution. The size of the deuteron polarizability corrections, including the 
elastic parts to Fig. is found to be of order 2 meV, and needs to be calculated to about 5% accuracy 
to allow a measurement to draw useful conclusions. 


There are a number of available polarizability calculations 14,24-28 that are mainly non-relativisitic 
and depend on potential model considerations of the deuteron bound state. The most recent includes 
the elastic contributions and quotes an uncertainty limit well within the requirement. In addition, much 
of that calculation is supported by a calculation 28 that uses the zero-range approximation which 


indicates that much of the result is relatively model independent. 
The result from 26 is 


AEJp^s = 1.665(16) meV , 


(26) 


where we have updated a small term coming from hadronic corrections from excitations of the proton 
and neutron, and left the uncertainty limit as it is found in (2^ . 

One would like in addition a calculation by a purely relativistic method analogous to the calculations 
done in 15 -|^ for the proton case to verify the results. These calculations involved using data on 
electron proton scattering plus dispersion relations to obtain the Compton amplitudes that make up 
the lower part of the box in Fig. Changing to the deuteron case requires data on inelastic and or 
quasi-elastic deuteron breakup which is not well enough known in the decisive kinematic range. The 


calculation has nonetheless been attempted 29 , with the result coming with a large uncertainty limit. 


AE 2 p _25 = 2.01(74) meV. 


(27) 


The better news here is that the same experiment that will better measure the deuteron charge radius 
using electron scattering will also gather useful quasi-elastic data that could help reduce the uncertainty 


in Eq. (27) to a useful level 23 
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Figure 6: One photon exchange. 


4 Calculating the proton radius dependence 


A number of theoretical calculations and corrections are needed to convert the measured energy split¬ 
tings into results for the charge radius. The original publication, for example, gave a list of more than 
20 corrections, albeit only 3 of them were large enough to affect the hnal result noticeable; Ref. 
also contains tabulations of the calculated theoretical corrections. 

There also other corrections. A year ago one might have worried that some of the corrections 
were single-sourced, i.e., had been calculated by only one person or group. Today, any thinkably 
important correction has been multiply checked; some reviews and/or additional small corrections 

We shall mention here only a few topics, mainly the higher order two- 


31, 32 


are given m 

photon corrections to the spectroscopy (because of its relevance to the discussion of the sources of the 
discrepancy), the two-photon corrections to scattering measurements (because they are still somewhat 
unsettled, though not apparently a decisive problem), and the lattice gauge theory proton radius results 
(because they may be very relevant in the future). 


4.1 Two photon corrections to the energy levels 


The two photon corrections were hrst calculated non-relativistically by Friar in 1979 14 , whose result 
added to the lower order term was given as 


A-B = (fll - 


(28) 


where R^ 2 ) is the average separation-cubed of the proton charge distribution. 


-^^2) = / Sr 2 PeSi) |fi - r2| 


W(’^2) 


(29) 


and Friar has referred to R{ 2 ) as the third Zemach radius. 

Today the calculation should be done relativistically and in momentum space, as pioneered by 
Pachucki 33 . The higher order perturbation theory calculation is represented by the diagram in Fig. 
and Friar’s result can be shown to be the non-relativistic limit of the elastic part of Fig. meaning the 
part where the blob in the Figure is just a single proton. 

An updated re-evaluation 17 has both updated this contribution and quantihed the error limits on 
its numerical value, using the latest input for what is known about the blob in the Figure. Neither this 
evaluation nor other related ones obtain a result large compared to 320 peV [l6p8p4fp^ . More detail is 
pushed forward to Sec. 5.1, where it is needed to explain some controversy regarding a contrary claim (3^ 


that off-shell contributions existing within the blob in the Figure could indeed give contributions as large 
as the 320 /leV needed to resolve the proton radius problem. 
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4.2 Two photon and Coulomb corrections in scattering 


There are corrections to electron-proton scattering coming from mnltiple photon exchange. These are 
Coulomb corrections, and are dominated by two photon exchange. The corrections at very low energies 


and very low momentum transfers were worked out long ago 38[ 39 , and at higher energies and higher 
momentum transfers in more modern times 


40 43 


which 
but the 


The experimenters have made two photon corrections, basically using the older works 38,39 
used pointlike protons. There has been some discussion of details of the corrections 44-46 , 
discussions did not center on the assuming of pointlike protons in the corrections. 

Recently, Lee and Milstein (see also 48-^) have noted that two photon corrections to the 


proton charge radius themselves have a correction due to the proton charge radius, and while not 
claiming this extra correction is large enough to explain the proton radius puzzle, they do claim it 
is too large to neglect. Including the extra correction reduces the proton charge radius measured in 
electron scattering by about 0.8%, for the lowest Mainz electron energy of 180 MeV and taking the 
—)■ 0 limit. The paper 47 exchanges only Coulomb photons, and it is unclear if the form factor 


they insert is Dirac or electric, but the basic statements appear correct. More recently, Gorchtein 
has included some relativistic corrections to the Lee-Milstein result and also given a possibly useful 
expression for the low-Q^ contribution from inelastic intermediate states in two-photon exchange. 


4.3 Proton radius ab initio calculations 


Ah initio in the context of QCD includes lattice gauge theory calculations. 

The current state of lattice calculations for nucleon radii and nucleon form factors is limited because 
lattice simulations with physical quark masses are computer time intensive and because lattice com¬ 
putations with disconnected quark lines are also time intensive. One escapes the hrst by calculating 
with heavier quarks, indicating the quark mass by giving the pion mass, and escapes the second by 
only studying isovector nucleon quantities. Isoscalar form factors do require disconnected diagrams, 
diagrams with quark loops not connected to the quark lines emanating from or ending on the lattice 
nucleon source or sink. The gluons that attach the quark loops to the valence quarks are not indicated 
in lattice diagrams, hence the phrase “disconnected.” However, the disconnected diagrams contribute 
equally to proton and neutron, and the isovector case has been assayed. 


One may specihcally focus on nucleon radii calculated from lattice gauge theory. In the future, it 
may be possible and desirable to calculate using a dedicated correlator which gives directly the slope 
of the form factor at zero momentum transfer. Finding such correlators by taking derivatives of known 
correlators is studied in 52 for lattice calculations of form factors that multiply Lorentz factors that 
go to zero at interesting points. Applications in (5^ are to form factors for semi-leptonic scalar meson 


decay, and to hadronic vacuum polarization corrections to the muon (^f — 2). 

At present, lattice calculations of nucleon radii proceed by calculating the form factor at several 
non-zero htting to a suitable form, typically a dipole form, and hnding the radius by extrapolating 
to zero Q^. Results are available only for the isovector nucleon 53 -55 . Ref. 54 presents a plot of Dirac 


radius results for lattice calculations at various pion masses. The results are typically about 50% of 
the experimental (electron measured) isovector Dirac radii-squared, albeit rising with decreasing pion 
mass, with quoted uncertainties in the 10% range. One may say there opportunity for further work. An 
uncertainty of 1% or less for the proton alone is needed for a lattice calculation to impact the proton 
radius puzzle. 
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5 Explanations 

If the muonic Lamb shift proton radius is correct, three explanations have been discussed to explain 
the discrepancy: 

1. There are unexpected QCD corrections that can give ~ 320 /reV energy to the muonic hydrogen 
Lamb shift. 

2. There is some “new physics,” a so far undiscovered new particle or particles with new interactions 
that violate e-/x universality. 

3. There is a problem with the radius determined from the full collection of electronic measurements. 
We will discuss each of them in turn. Each of them may work, but each of them has problems that 

we shall try to mention clearly. 


5.1 “Haywire” hadronic behavior 


It is possible that unexpected QCD effects could exist with unexpected magnitude. The context of the 
discussion is the two-photon corrections that have already been diagramed in Fig. The generally 
accepted value has already been quoted, Eq. (20), and the standard analyses are in general agreement 


on this value 15 18,34 . However, there are parts of the correction that are strictly speaking unknown. 


The standard estimates of these parts have them extremely small, but perhaps it is otherwise 56 and 
an examination of the issues will be given here. 

We shall give some detail in order to explain how the discussion arises. Relativistically, the cor¬ 
rections from two-photon exchange are obtained from the Feynman diagram Fig. On the scale of 
nuclear effects, the Fermi momenta of the atomic electrons or muons can be neglected even in a precision 
calculation. The bottom part of the diagram is forward, off-shell Compton scattering. The intermediate 
state in the lower leg may be a proton (the “elastic” case) or may be more than a proton and then 
“inelastic”. The lower part of Fig. both elastic and inelastic, is described by the Compton tensor, 
which is 




SttM 




= 1 ^ ] TM Q^) + ^ (P" - 


P 


p-q 


g" T2 (z/,q2), 


(30) 

g^ y ' ■ \ j \ g^ y 

where = —g^ and v = p ■ q/M (Fig. |^. Straightforwardly, the energy shift in terms of Ti and T 2 is, 

f + ^Ql)Ti{tQo,Q^) - {Q^ - Ql)T2{tQo,Q^) 

/ a ty n4('n4 I /I A/f2n2'i 


AE = 


TT 


(31) 


where we have done a Wick rotation so that go = iQo and Q = q. 

The serious problem that we don’t know what Ti and T 2 are. But we do have the optical theorem, 
that the imaginary part of a forward scattering amplitude is related to the total (elastic plus inelastic) 
cross section for the same incoming states. 

In this case, the total cross section is the the one for electron-proton scattering. Fig. The cross 
section over a range of kinematics has been measured at SLAC, DESY, Bonn, JLab, Mainz, and other 
laboratories and codihed in the spin averaged case in terms of two structure functions Fi(z/, Q^) and 
F 2 {i'^ Q^). These then give the imaginary parts of the Compton amplitudes as 


ImTi(z/,QQ = 


( 32 ) 
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Figure 7: Electron-proton scattering. 


(the normalization of the Fi is standard, the normalization of the Ti follows Eq. (30)). From the 


imaginary part of the Ti and the Cauchy integral formula, leading to what are known as dispersion 
relations, we can obtain the full amplitudes. 

But there is a glitch, and this glitch underlies the present discussion. When using the Cauchy 
integral formula, we are integrating in u and some of the integrals do not converge at high u. This 
affects Ti but not T 2 . So for T 2 we can write a simple dispersion relation, but for Ti we instead write 
a dispersion relation for Ti/u. The convergence is now fine, but now we have a pole a.t u = 0, and the 
price is that we need to know Ti(0,Q^). The dispersion relations are 


and 


Ti(go,Q')=Ti(0,g2) + 


T2{qo,Q^) = ^ 




2'kM 








^0 




<& 


(33) 


(34) 


The dispersion relation for Ti is a subtracted dispersion relation and the term Ti(0,g^) is the “sub¬ 
traction function.” 

To get the energy shift, there is still an integral needed over . Hence we need to know Ti(0, Q^') 
for all non-negative . This requires modeling because only the point = 0 corresponds to physical 
kinematics. The slope in of T(0, Q‘^) a.t = 0 can be obtained from low energy expansions 
and more recently the curvature at = 0 has been obtained using chiral perturbation theory 


Additionally, the scaling behavior at high-Q"^ is hxed from perturbative QCD considerations 57 . Hence 


17 


18 


it is possible to make a reasonable, to most workers, treatment of the subtraction function and assign 
it a reliable uncertainty limit. Inserting inelastic and elastic experimental knowledge of Fi and F 2 
for other terms in the integrals is straightforward. Using up to date data and treatment of the 
subtraction function from 18 one has (—4.2(1.2)-|-12.7(0.5)-|-29.5(1.3)) /reV from the subtraction 


term, the inelastic contributions from Fi^ 2 , and the elastic contributions from Fi^ 2 , respectively, summing 


to the result quoted earlier in Eq. (20). The result is noticeable (and, incidentally, gives a large fraction 
of the total uncertainty quoted for the muonic Lamb shift measurement), but it is not the 320/reV 
needed to explain the discrepancy while keeping the CODATA proton radius 

However, it is still true that the subtraction function is not really known. It is possible to speculate 
on forms for Ti(0,g^) that are innocuous for low but have a large peak at (say) > 10 GeV^ 
before settling down to the correct asymptotic behavior. It is possible to thereby obtain 320 /reV for 
the subtraction term, and so “solve” the proton radius problem as has been done in 56 . 

But the modihcation of the subtraction function does have further consequences. One should think 
about the Cottingham formula 58 -60 , which gives the electromagnetic contribution to the neutron- 
proton mass difference. The Compton tensor also enters here, illustrated in Fig. where the two 
photon lines that attach to the nucleon are tied together rather than attached to a lepton. One has 
the same subtraction function. The published subtraction function proposed to solve the proton radius 






















Figure 8: Electromagnetic contribution to nucleon mass. 


problem contributes a startling 600 MeV to the individual nucleon terms in the Cottingham formula. 
The measured neutron-proton mass difference is 1.29 Mev 10 . Even if one makes the large term in the 


subtraction function isoscalar, so that is does not affect the mass difference, one is still at this moment 
left with an implausible 600 MeV electromagnetic contribution to the individual nucleon masses. 


The validity of perturbative QED for electron-proton interactions has also been questioned 61,62 


but these explanations would seem to require more explanation of how the validity of perturbative QED 
can be vitiated in interactions with hadronic bound states. 


5.2 Exotic Explanations 

A possibility is that the discrepancy is a signal of new physics, a new phenomenon, perhaps an exchange 
of a new particle, that we are uncovering for the hrst time. See 63 -72 . Some of the references discuss 


new physics possibilities, some discuss problems that ensue, some discuss both, and some suggest 
possible further experiments. 

shows that a new attractive (the 320 /reV shift being downwards) particle that interacts 


Ref. 73 


equally with electrons and muons would in fact lead to a larger radius from muonic hydrogen, the 
opposite of what is observed. Hence we need to consider a breakdown of muon-electron universality. 
A difference between the interactions of the muon and electron, other than the obvious gravitational 


one, has been long sought and not found (some early hints 74,75 were considered “not experimentally 
signihcant” even in the original sources). The overall negative result of these searches will lead to 
difficulty in the present context. 

In particular, enhanced coupling to muons necessarily affects the theory for {g — 2)^, where there is 
in fact a known discrepancy, albeit one that is small in percentage terms compared to the proton radius 
discrepancy. However, without hne tunings, the effect of any new particle that can explain the proton 
radius puzzle will have a large effect upon (^f —2)^. Hence there is a sticking coupling: one cannot claim 
the theory corrections to (^f — 2)^ are under good control unless the proton radius puzzle is understood. 

Considering the workings of a new particle in a way that is common to (at least) Refs, 
consider a particle that couples to muon and protons, but with no coupling or weaker coupling to 
electrons. The particlecan be scalar or vector. Pseudoscalars or axial vectors give small contributions 
at low momentum transfers, i.e., to atomic binding, for similar couplings. 

For illustration, the scalar case has an interaction Lagrangian 


Cs = -C"^ it> II - cl p 4) p, 


(36) 


where muon and proton helds are denoted by the symbol for the particle and 0 is a scalar held, and 
cr are constant coupling parameters. This gives a non-relativistically calculated 25* -2P energy shift 


AE= r‘^dr V{r) {R 
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Figure 9: New particle scalar couplings vs. mass. 




Figure 10: Corrections to the muon magnetic moment due to new particle exchange 


where M is the mass of the scalar 0 and = m^mp/{m^ + rrip). The (j) mass hxes the product of 
the couplings, since we must obtain AE = 320 /reV. For the case Cg = Cg = Cs, the result is shown 
graphically in Fig. A plot for the vector case would look the same. 

Now for = {g —2)^/2, the discrepancy between theory and experiment is 2 part-per-million (ppm) 
using only standard physics. In more detail, the results for are 


a^(data) = (116 592089 ± 63) x 10“^^ [0.5 ppm], 

a/,(thy.) = (116591840 ±59) x 10“^^ [0.5 ppm], 

= (249 ± 87) X 10"“ [2.1 ppm ± 0.7 ppm]. (37) 


(As a note, the above is from 76 , where the completed 5-loop QED {g — 2) calculations were reported. 
Their theory number used the hadronic corrections from Hagiwara et al. [^, which was one of three 
papers reporting hadronic corrections to {g — 2)^ in 2011. Hagiwara et al. was the largest, and in units 
of 10“^^, Davier et al. 78 were 20 units smaller, and Jegerlehner and Szafron 79 were about 20 units 
smaller still. Hence the discrepancy could be about 40 units larger than shown above.) 

What makes it possible to have small corrections to [g — 2)^ is the fact that corrections to the 
anomalous magnetic moment. Fig. IT, from scalar and pseudoscalar particles have opposite signs. The 
same is true for corrections from polar vector and axial vector exchanges. 

If, as an example, we include both new scalar and new pseudoscalar particles, we have the amended 
Lagrangian, 

Cs = -C^ (psfil^-C^4>sPp- (t>p - iCp 0 0p750p • (38) 

For a given mass and the already determined Cs, one can find the Cp required to keep at or below 


the known discrepancy. The result is shown in the left hand panel of Fig. 11 Another way of expressing 
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Figure 11: Top left panel: New physics S', P couplings required to both explain the proton 
radius problem and hne tune [g — 2)^. The solid line is Cg and the dashed line is Cp. Top 
right panel: Degree of hne tuning required for the scalar-pseudoscalar case. Bottom panel: 
New physics V, A couplings required to both explain the proton radius problem and hne 
tune (g — 2)^. The solid hne is Cy and Ca is dashed. 


the result is in the degree of hne tuning required. The diherence between the individual contributions 
to the anomalous magnetic moment from the scalar and pseudoscalar exchanges, compared to the 
contribution from either one of them, is shown in the center panel of Fig. 11 The hne tuning becomes 


extreme as the mass of the exchanged particle rises. However, at low masses the contribution of the 
scalar by itself to {g — 2)^ is small enough that no cancellation from a pseudoscalar is needed. This in 
fact is the situation envisioned in 63 . A plot showing the couplings needed for the vector-axial vector 


case is shown in the right-hand panel. Fine tuning is also required here, but not to the degree seen in 
the scalar-pseudoscalar case. 


5.2.1 Constraints from iF-decay 

New particles that couple to muons can give radiative corrections to any decay involving muons, and 
an immediate example is radiative corrections to —)• gAu, Fi g, Depending on the coupling to 

electrons, the 0 particle in Fig. 12 can either exit the detector 65 69[|80|, or decay into e’''e“ pairs 


81 


Both possibilities are interesting. 
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Figure 12: Kaon decay with an extra neutral scalar, 0, either scalar or pseudoscalar. 


If the (j) leaves the detector without decay or interaction, the only visible decay particle will be a 
muon, with energy below the energy it would have in the two body decay —)■ Experiments 

have hunted for —>■ /i^ + invisible 80 , where “invisible” means particles that are neutral but not 


photons. The original motivation appears to have to explore the possibility that there were multiple 
neutrinos in the final state, but the bound applies to the present case also. The results lead to mass 
constraints that eliminate low mass (m^ < lOOMeK) new vector particles or eliminate middle mass 
(90 < < 200 MeV) new scalars 

If the 0 can decay into e+e“, the bounds of the last paragraph don’t apply. The decay that 
does occur, K —)■ /iz/0 —)■ competes with the standard QED radiative decay K —/rz/y* —)■ 

The former case of course has the mass concentrated at the mass of the 0 and with 


82,83 


/iz/e^e 

a mass choice, the couplings to the muon are hxed. The 0 is certain to decay in the detector unless the 
electron couplings are very small, and theory studies show a dehnite peak above background 81 
if the new particle in fact exists. 


5.2.2 HFS and other constraints 


Ref. points out that the new particle exchange can contribute to the HFS in muoninum [the fie 
bound state], which is accurately measured and for which the standard model theory is fairly accurate . 
Hence a small new contribution can matter, and and new particle that couples to muons can still couple 
to electrons through a muon loop that attaches to photons that carry the interaction to the electrons. 
The size of the contribution requires the new particle mass to lie under 25 MeV. 

Further, the measurement of the HFS in muonic hydrogen [^, 

AEnfg = 22.8089(51) meV, (39) 


agrees within uncertainties with the value calculated using known physics 20 . This also places pressure 


on exotic proton radius puzzle explanations. Pseudoscalar or axial particles were introduced for the fine 
tuning of {g — 2)^ and these particles do not contribute strongly to the overall muonic hydrogen energy 
shift. But they do have in the nonrelativistic limit a noticeable spin-spin force and hence affect the 
HFS predictions. We will give a few details. To state the conclusion at the outset, lower mass exotic 
exchanges are still possible, high mass exotic exchanged particles are not. 

For definiteness, consider the axial exchange in Fig. 
and Sullivan 


13 The calculation is similar to one of Drell 


84 in another context. In the nonrelativistic limit the matrix element is. 


Ma = 2m^ 2mp \ 


1 


and the contribution to the HFS for the 2S state is 

2mla^ m\{m\ + 




dvr 


m. 


(rriA + m^a)^ 


(40) 


(41) 
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Figure 13: Axial vector exchange in an atomic system. 




Figure 14: Contributions of exotic particle exchanges (in one model [^) where the exotic 
particle is axial vector (left panel) or pseudoscalar (right panel). To avid conflict with data 
and with standard model calculations, the contribution should be below 5.1 /reV, indicated 
by a the dotted line. 


With as a function of the exchanged mass already obtained from the fine tuning constraints, 

there follows the HFS results vs. mass shown in Fig. 14, left panel. To obtain a contribution to the 
HFS that is below the allowed tolerance, the axial vector mass should be below about 13 MeV. (The 
corresponding result for the polar vector exchange is also shown on the plot.) This, of course, has been 
done in one particular model 85 , but other exotic models explaining the proton radius puzzle will 


behave qualitatively similarly. 

A similar plot for the pseudoscalar case appears in Fig. 14, right panel. In this case the allowed 
mass can be somewhat higher, but still needs to be below about 35 MeV. 

Another series of constraints comes from searches for decays T —)■ 7 + invisible or J/'ip —)■ 7 + invisible, 
which will emphasize the targeted couplings of any putative new particles. The first of these is a bb 
bound state, and the other a cc bound state. The motivation for these searches was to find glueballs. A 
large part of the decay at the perturbative level should be into jgg final states, where is a gluon, and 
it was hoped that the gg would often bind into glueball state with a definite mass an hence bumps in 
the 7 spectrum. No such states were found. The data left behind is a serious constraint here, because 
the couplings for a given mass are now known for the proton and if the coupling were the same for the 
5 or c quark, the new particles that £x the proton radius problem would be seen here 
leads to forbidding coupling of the new particle to second or third generation quarks. 

Finally, Ref. 


64 . This then 


72 


emphasizes the need for a gauge-invariant, renormalizable, theory for the new 
physics explanations. The interactions that have been written should be viewed as low energy ap¬ 
proximations to a complete theory. A non-renormalizable theory has interactions whose effect grows 
excessively with energy, and in the present case the high energy problems are visible in the large radiative 
















corrections toW^^v decay 72 . 


Alltogether, the prospect of new particles explaining the proton radius problem is difficult but not 
impossible. Low masses are favored. One can to “explain” two vexing problems of modern physics, 
the proton radius problem and the muon anomalous magnetic moment problem, using the existence 
of several possible particles with unusual properties. The price is high. There need to be targeted 
couplings and hne tuning. Nonetheless, the idea is not eliminated. 


5.3 Electronic scattering measurements 

There are many experiments using electrons that determine the proton charge radius, and only one 
using muons. So it is natural to focus on the muon experiment, but one could also entertain the idea 
that the charge radius determined using electrons is not as secure as the error limits suggest. We begin 
by discussing the radius obtained from electron scattering and list some alternating analyses. 

One has begin by respecting the £t of the experimenters themselves to their own data . 

The quality and number of electron scattering data from the 2010 Mainz experiment is remarkable. 
They have 1422 data points, with from about 0.004 GeV^ to about 1 GeV^. They £t their entire 
data set with a number of different parameterizations of the electric and magnetic form factors, and 
this is where one obtains the electric radius quoted earlier. 

However, the dispersive analysis reported in 88 also uses only electron scattering data but obtains 
a charge radius in accord with the /r-H value. They treat the form factors as analytic functions, and 
there are cuts that show that hts to the data using functions like polynomials in are not guaranteed 
convergent for above 4m^ (on this point see also 89,^). They extend the analysis to include 


timelike data, which can done using the analyticity properties of the scattering amplitudes and the 
dispersion theory. One viewpoint on the dispersive treatment is that it gives £t functions that satisfy 
the necessary analyticity requirements at all an hence are convergent. The results of 88 give a 
proton charge radius Re = 0.84(1) fm. 

On the other hand, the authors of 89 also respect the analytic properties of the form factors, but 
they hnd a result well in accord with the GODATA value. They conformally map onto another 
variable, where the Taylor expansion in the other variable is expected to be convergent everywhere. 
They £t the data catalogued in |^, which clearly does not include the newest data, and find stable 
fits with Re = 0.870(26) fm (combining in quadrature their two quoted uncertainties). An update of 
this work including all currently available data is underway 91 . However, currently available is similar 
work using the conformally mapped variable done by 


90 


htting the entire Mainz data set and hnding 
the lower Re = 0.84(1) fm proton radius. A difference between 89 and 90 , in addition to the different 


data sets, is that the former restricted the size of the coefficients in the expansion, whereas the latter 
did not. The former may argue that the coefficients are thereby a physically more reasonable size, the 
latter may point out the signihcantly lower per degree of freedom obtained. 

Going in another direction, the charge radius is a low quantity, so that low data can suffice to 
determine it. The determination requires hnding the difference of the form factor from unity at small 
but the measurement range should not be too small lest uncertainties in the measurement make 
determination of the slope impossible. On the other hand, for a data selection covering a wide range 
of the form factors measured at the high part of the range have little relevance to determining the 
slope at = 0. Ref. 92 suggests that the range of needed to determine the slope is in the range 


0.01 to 0.06 GeV^. 

Gonsidering this hrst over a somewhat narrower range, of the 1422 Mainz 2010 data points, about 
200 are at below 0.02 GeV^, and have small uncertainty. Fitting just these circa 200 data points 
leads strikingly to a radius agreeing with the /i-H value with about 2% uncertainty [Re = 0.842 (20) fm]. 
An additional detail is that the £t has terms up to (Q^)^, terms beyond that being not determinable 
from and by expectation not important in low data. The curvature, the coefficient of the term. 


20 













The attitude one takes to this depends on 
The result that the radius comes out to the lower 


obtained has small central value but large uncertainty limit 
whether one focuses on the former or the latter 
value continues even as the upper limit of in the data selection doubles or triples 
Also hnding large radii are Ref. 92 


93 


points out that the proton matter density at large 
r might give an effect on the form factor that is small overall, but large at very low and vitiate the 
extrapolation to = 0. (An similar point was pursued in 95 , but with an implementation giving 
form factors far from data at measured 


22,^.) They also, perhaps not in concert with the first 


point, advocate including high data when fitting the charge radius. In any case, fits are done with 
continued fractions 11 and sum of gaussians forms (including the Mainz data), with a result in 


the more recent work Re = 0.886(8) fm. 

Perhaps the only sure conclusion from these conflicting claims is that the uncertainty in the proton 
radius obtained from electron scattering is larger than what is quoted. The large number of recent 
analyses on related issues illustrates that the precise htting of the proton form factor in the low 
region is not straightforward and remains a matter of ongoing discussion. Data at still lower and 
muon scattering data would be helpful, and there are such plans for both e-p and p-p elastic scattering, 
which will be mentioned in separate sections below. 

Of course, electronic measurements of the proton radius come from atomic spectroscopy as well as 
from electron scattering. Again, no single measurement gives the quoted better than 1% accuracy, but 
rather there are on the order of 15 measurements with accuracies of a few percent that legitimately 
give the smaller error limit when combined. If one wants to entertain the idea that the proton radius 
is really the lower muonic Lamb shift value, one will have to explain why the spectroscopy results give 
the larger value. However, there are no known problems with the atomic spectroscopy measurements 
and no alternative analyses obtaining different radii. 


6 The Future 


The proton radius puzzle remains unsolved. Further, it has implications for other fields. Since exotic 
explanations must impact other purely muonic processes, the theory for the {g — 2)fj_ measurements 
cannot be deemed under control until the proton radius puzzle is solved 


72 


There are a number of experiments relevant to the proton radius puzzle, and a small and perhaps 
incomplete catalog is given here, along with some rather brief comments about what the experiment is 
and why it connects to this puzzle. It is exciting to anticipate the large amount of potentially useful 
information that is coming. 


6.1 New CREMA measurements 

The published CREMA measurements are for muonic hydrogen [^[^. The collaboration also measured 
2S-2P transitions in muonic deuterium simultaneously with the first muonic hydrogen run, which are 
presently being analyzed, and a similar measurement of the Lamb shift in muonic helium ions has 
been performed more recently at the PSI. A combined analysis of these experiments will help to clarify 
whether the problem arises in the muonic sector, or if it is related to the proton, or the Rydberg 
constant, or if it originates from bound-state QED, or if it is something still more startling. 

6.2 MUSE 

The MUon proton Scattering Experiment (MUSE) at the Paul Scherrer Institute is a simultaneous 
measurement of p'^p and e^p elastic scattering, as well as of p~p and e~p elastic scattering. The main 
interest here is in filling in a category of measurement. That is, the proton radius has been measured 
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Figure 15: Initial state radiation in elastic electron-proton scattering. 


consistently using electrons in both scattering and spectroscopy; it has been measured using muons 
only via spectroscopy. This experiment will be the hrst see if muon-proton scattering gives the same 
answer as muonic hydrogen spectroscopy. 

In addition, the differences between positive and negative charge lepton scattering measure two- 
photon exchange effects, which are higher-order corrections to the scattering process. The experimenters 
plan to determine the proton radius from fip scattering at the level of about 0.01 fm, similar to previous 
e-p measurements. They hope to measure momentum transfers in the range 0.002 to 0.07 GeV^. There 
have been some test runs already, with two 6-month production runs scheduled for 2016 and 2017. 


6.3 PRad 

This experiment will measure e-p elastic cross sections at high precision and at very low four-momentum 
transfer squared, Q^, from about 2 x 10“^ to 2 x 10“^ GeV^, in Hall B at Jefferson Lab 98 . They 
obtain the low by having low scattering angles, which is possible because they have no bulky magnetic 
spectrometers, but just a target and a downstream lead-glass wall acting as a calorimeter. The absolute 
value of the e-p cross sections is calibrated by a well known QED process, Mpller scattering or e~-e~ 
elastic scattering, which will be continuously measured in this experiment with similar kinematics and 
the same apparatus. With accurate data and low Q^, they hope to obtain a sub-percent and essentially 
model independent extraction of the proton charge radius. 

6.4 ISR form factor measurements 

An experiment at MAMI is underway, aimed at measuring proton form-factors at very low momentum 


transfers by using a technique based on initial state radiation (ISR) [^. The idea is that ISR, Fig. 15 
degrades the energy of the incoming electron so that the momentum transfer to the proton can be 
quite low. The outgoing electron angle and energy are measured as usual, and with some modeling the 
diagram shown can be isolated from other radiative corrections and an accurate form factor obtained 
with as low as 10““^ GeV^. There was a pilot measurement in 2010, the full experiment has run in 
2013, and the analysis is continuing. 


6.5 Electron deuteron scattering 

Electron-deuteron scattering data currently available does not by itself give a accurate enough deter¬ 
mination of the deuteron charge radius to assess any discrepancy occasioned by the muonic deuterium 
Lamb shift. Fortunately, the hydrogen-deuterium isotope shift allows using the proton radius measure¬ 
ment to obtain the deuteron radius to high accuracy 100 , as discussed in Sec. 3.1 above. However, 


conhrmation from direct scattering measurement is underway at Mainz 
the analysis is ongoing. 


23 . Data has been taken, and 
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6.6 High precision Lamb shift in e-H 

This is a proposed measurement of the ordinary atomic hydrogen n = 2 Lamb shift using the Ramsey 


method of separated oscillatory helds 101 . This new measurement, with an anticipated uncertainty 5 


times more accurate than existing measurements, along with existing precise atomic theory calculations, 
will by itself allow for a new determination of the proton charge radius to an accuracy of 0.6 percent. 
Currently the quoted uncertainty in the proton charge radius from atomic spectroscopy comes from 
averaging over a large basket of measurements. Here one will get a lower uncertainty result from just 
one measurement. 


6.7 New measurements of larger e-H splittings at Garching 

Work is underway at the Max Plank Institute for Quantum Optics in Garching on a new precision 
measurement of the 2S'-4P transition in atomic hydrogen. This measurement, paired with the very 
accurate 15-25' measurement can give very accurately both the Rydberg and the proton radius 102[|103 


Currently, as relevant also above, the individual electronic hydrogen spectroscopy measurements or pairs 
of measurements have proton radius uncertainties of more than 2% and the small overall uncertainty 
comes from averaging over many measurements. Here one will get and low uncertainty results from one 
measurement. 


6.8 New measurements of larger e-H splittings at LKB, Paris 


Also in the category of new measurements of larger e-H splitting is work underway on a precision 
measurement of the 15-35 splitting in atomic hydrogen at the Laboratoire Kastler Brossel (LKB) in 
Paris. Again, this measurement by itself paired with the very accurate 15-25 measurement will give 
very accurately both the Rydberg and the proton radius. To get an idea of the precision necessary, the 
15-35 splitting is about 2.9 x 10^^ kHZ, and the difference between the predicted splittings with the 
CODATA radius and the p-iL radius is about 7.2 kHz or a relative difference of about 2.5 parts per 
trillion (ppt). A 1% accuracy for the proton radius will require a measurement to 1.8 kHz (0.6 ppt) or 
less. The 2010 result from the LKB for the same splitting, which will be superseded, had a uncertainty 


limit of 13 kHz 104 


6.9 Alternative measurements of the Rydberg 

As discussed earlier, the Rydberg constant and the proton charge radius are determined in the same 
spectroscopy experiments. If the muonic hydrogen Lamb shift proton radius is used, the Rydberg 
constant would change by 4 standard deviations. It may be possible to independently measure the 
Rydberg constant by using Rydberg states (very high-lying states) where the effect of the proton size 
is negligible. Consideration of this possibility is underway at the National Insititute of Science and 
Technology (NIST, USA) [TOS] . 


6.10 Trumuonium at JLab 

True muonium is the ii~ bound state. If there is a exotic exchanged particle that couples to muons, 
it will affect the binding energy of true muonium. The Heavy Photon Search experiment (HPS) at Hall 
B at Jefferson Lab will as a side project also have the possibility of finding true muonium |106| . There 
is hope of seeing 60-100 true muonium events; it is of course some distance from just observing true 
muonium to measuring its spectrum accurately. 
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7 Conclusion 


The proton charge radius puzzle remains a major problem. 

The conflict is between measurements using electrons and the one muonic experiment, which mea¬ 
sures the muonic hydrogen Lamb shift. On the other hand, the muonic experiment has by far the 
smallest uncertainty limit of any proton radius measurement and while it is a difficult experiment to 
make work, it is does not seem to be a difficult experiment to analyze. Further, the theoretical correc¬ 
tions have been reexamined and repeated by a number of worker since the announcement of the result, 
with no generally accepted problem has been found. 

Perhaps the discrepancy is due to new physics. Perhaps the explanation is an ordinary physics effect 
that has been missed. Perhaps, the muonically measured radius will come to be the accepted number. 
Time will tell. 
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